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I. INTRODUCTION 

One cannot overstate the practical importance of 
understanding quantum decoherence. The remarkable 
emergence of quantum engineering and the pursuit of 
quantum computation have resulted in the creation of 
entire industries relying upon the minimization of the ef- 
fects of decoherence. 

Less practical but perhaps more significant is the in- 
sight to be gained into the theory of fundamental pro- 
cesses. Many have sought a natural explanation for this 
experimental signature of an arrow of time for quantum 
systems, which have states represented theoretically by 
density operators with intrinsically reversible evolution 
in time [341 ■ 

Phenomenologically, one imagines irreversible behav- 
ior to result from uncontrolled interactions between an 
experimental system and its environment. An arrow of 
time is thus usually accommodated in quantum theory 
by modeling the system of interest as a subsystem em- 
bedded within a quantum mechanical environment Q . 
While the time evolution of the environment plus subsys- 
tem is unitary, the time evolution of the subsystem alone 
need not be. In practice, one typically neglects memory 
effects of the environment d, 01 and models this reduced 
dynamics of the subsystem using a Lindblad-Kossakowski 
master equation. 

Nature has provided us also with decaying states and 
scattering resonances, which are thought to be signatures 
of an asymmetry in time Q even for closed systems. This 
time asymmetry is called intrinsic to distinguish it from 
irreversibility and the extrinsic arrow of time, which is 
consequent to environmental interactions. It is endowed 
at the microphysical level and is thus understood to have 
a different physical origin [9]. If one demands a rigor- 
ous, theoretical unification of scattering resonances with 
decaying states, then one requires time evolution gener- 
ated by a semigroup [l0|. Semigroup time evolution is a 
theoretical expression of microphysical time asymmetry. 

Kossakowski and ReboUedo have recently begun to ap- 
ply concepts from the intrinsically time asymmetric the- 



ory to the formulation of a new master equation 
Their approach is based on the similarities between the 
energy spectra of subsystems embedded in an environ- 
ment and of scattering resonances embedded in a contin- 
uum. 

Here we take a different approach. There is a direct re- 
lationship between the intrinsic time asymmetry and the 
extrinsic arrow of time. If one chooses microphysically 
asymmetric time evolution, one must very carefully con- 
struct the theoretical images of quantum systems. One 
also finds a new correspondence between time evolution 
parameters and the passage of time in the physical uni- 
verse. It is just this new correspondence, and not some 
mathematical effect, that leads automatically to predic- 
tions of decoherence for open systems. 

After explaining briefly in Section |lT] what we mean by 
intrinsically asymmetric time evolution, we will demon- 
strate in Section inil how one must represent experimental 
systems. 

In Section IIVI is an application. We derive a predic- 
tive probability for the measured decoherence of quan- 
tum systems undergoing Rabi oscillations. Furthermore, 
we match a general experimental result that has been 
particularly puzzling. Using our new approach, we con- 
clude that what has been called Excitation Induced De- 
phasing [T^l is a consequence of the indistinguishability 
of separate, uncontrolled interactions between quantum 
systems and their environment. 



II. GROUPS AND SEMIGROUPS 

The theoretical image of a quantum mechanical system 
is an operator algebra defined in a linear scalar-product 
space, $ Ts]. The vectors 4>n span the space $, and every 
linear combination of the 0„ can represent the state of 
the physical system. For simplicity, in this section let us 
consider a single vector, 0, that satisfies the Schrodinger 
equation. The vector </> spans a one-dimensional subspace 
of $, and this ray represents the quantum system in a 
pure state. 
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For the scalar-product space, one has historicahy 
chosen the Hilbert space: $ = "H. With this choice, 
Stone and von Neumann showed [3, [l^ that the time 
evolution of solutions is given generally by 

(j)(t) ^ e~ ^ 4> , -(X)<t<oo. (1) 

To require € H is to constrain 0, just as one con- 
strains the general solutions of any differential equation 
by enforcing proper behavior at the boundaries. For this 
reason, one calls € "H a choice of boundary conditions. 

The time evolution in (IT|) can be described by the one- 
parameter group of unitary operators 

U{t)^e-'^, ~oo<t<oo. (2) 

For every evolution, U{t), there exists the inverse, 
given by U{t)-'^ = U{-t). This means that 
U{—t)(j) is the time-reversed version of U{t)(l), and that 
both are solutions of the Schrddinger equation. The re- 
sulting dynamics is intrinsically symmetric in time. 

One often chooses state vectors from a different scalar- 
product space [IB], though usually without realizing it. 
When using Dirac kets {\E), \x), etc.) with continuous 
spectra of eigenvalues, one has chosen (p £ S C H, where 
the Schwartz space, S, is a subset of H [l3| . 

Though they are measured differently, one typically 
assumes that scattering resonances are physically equiv- 
alent to decaying states. To find a rigorous theory in 
which the two can be considered also mathematically 
equivalent, one must constrain the state vectors even fur- 
ther [13 : 

set of possible states {0} = $_ c H C (3) 

where $_ denotes the Hardy space of the lower complex 
semiplane, and <I>^ is its dual. Note that, in practice, ^ 
is not a limiting restriction. Any vector in the Hilbert 
space can be approximated with arbitrary precision by 
vectors in the Hardy space [ll]. This section contains 
only a sketch, and for details the interested reader is re- 
ferred to [1, [H, US] and and the numerous references 
therein. 

The only mathematical consequence of ^ we shall 
consider here is that solutions of the Schrodinger equa- 
tion become [ll] 

(j){t) = e"^0 , 0<t<oo. (4) 

Note the lower bound on the time parameter. This is no 
longer time evolution given by the unitary group of ©. 
The time evolution is instead given by the one-parameter 
semigroup of operators, 

f7(i) = e"^, 0<t<oo. (5) 

Being a semigroup means that the inverse, U{t)~^, of any 
element, U{t) with t > 0, does not exist. In other words, 
the time-reversed version ofU(t)(j) is no longer an avail- 
able solution of the Schrodinger equation. This dynamics 



is intrinsically asymmetric in time [9|. Prom a choice of 
boundary conditions, asymmetric time evolution is en- 
dowed at the microphysical level and is independent of 
any interaction of a system with its environment. 

The use of semigroups is not new to physics. First, 
expert readers will notice that ([S]) is not the semigroup 
often assumed for the family of dynamical maps for a 
reduced density matrix jSj. It has also been suggested 
in the past that one might describe relativistic, unsta- 
ble particles using representations of the Poincare semi- 
group [Hm. 

III. APPEARANCE OF DECOHERENCE 

With our choice of boundary conditions (jS]), and thus 
the microphysical time asymmetry, comes a handful of 
consequences affecting the representation of experimen- 
tal systems by their theoretical images. For open sys- 
tems, predictions of decoherence follow automatically 
from standard quantum mechanics. 



A. Coordinates and Parameters 

In a recent review of time in quantum mechanics [25} . 
Zch explains that, for the non- relativistic theory, one 
identifies the parameter t in the Schrodinger equation 
with "Newton's absolute," or coordinate, time. Because 
in the master equation formalism, one models experi- 
mental systems as subsystems of arbitrarily large envi- 
ronmental systems also evolving in time, consistency re- 
quires one to identify t as such. 

In some cases, however, it is known that such an iden- 
tification can be problematic. For instance, to formu- 
late a sensible version of the time-energy uncertainty re- 
lation, one must distinguish between "external time," 
which is measured by laboratory clocks, and "intrinsic 
time," which parametrizes the dynamical evolution of ex- 
perimental systems [26| . 

We will demonstrate that the correct application of 
the intrinsically time asymmetric theory requires one to 
distinguish between external and intrinsic time. What is 
called external time in [2^ we will represent with time 
coordinates. Time coordinates are the time part of the 
space-time coordinates. They are physically insignificant 
time labels for events, and their values are given by the 
laboratory clocks, which are not dynamically connected 
to experimental systems. We will label time coordinates 
with a tilde: t. 

What is called intrinsic time in [26j we must iden- 
tify with the t in the Schrodinger equation and its so- 
lutions (HI). We will represent this intrinsic time with 
time parameters, which parametrize the time evolu- 
tion of state vectors. Time parameters always correspond 
to durations, and we will continue to label them with the 
letter t. In Table U is a reminder of the distinction we 
will have to make. 
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Name 


Symbol 


Description 


coordinate 


i 


time labels for events; 
marked by clocks 


semigroup 
time parameter 


t 


parameter for Schrodinger 
equation and solutions (|4]) 



TABLE I: Time coordinates and parameters. 



Considering the time asymmetric boundary conditions, 
that one must distinguish carefully between coordinates 
and parameters is obvious because possible values of t 
and t are no longer chosen from the same interval: 

t€[0,oo) and ie(— 00,00). (6) 

In what follows, we will assume that one can infer from 
measurement or from a preparation procedure the na- 
ture of a density operator meant to represent the state of 
a physical system. Measurements and preparations are 
performed in the coordinate time, t, of the laboratory, 
so this inferred density operator is also a function of t. 
Because only time durations are significant, we will call 
it p(t — tprep), where tprep is the time coordinate value 
when the state of the system is prepared (see below.) 

Any density operator written as p{t) is understood to 
be a function of the time parameter, t, of theory. It is 
calculated using standard quantum mechanics. And as 
explained below, for a given system one cannot always 
equate p(t — tprep) with p{t). To summarize. 

Calculated Inferred from measurement 

pit) p{i- tprep) 



B. Preparation Time and t — 

In the theory, there is now a distinguished value, t — 0, 
of the parameter used in Q to parametrize the time evo- 
lution of state vectors. This distinguished value of time is 
a mathematical feature that is phenomenologically signif- 
icant and identifiable in comparisons with experimental 
data. 

One defines a physical observable by a prescription for 
how it is to be measured. In the theory, one re pres ents 
observables with linear operators in a space $ [l3|. In 
the Schrodinger picture, the Born probability to find the 
observable A = in the state p(t) — |0(i))((/)(t)| is 

VA{p{t)) = Tr(Ap(t)) = m<f>m^ 0<t<^. 

(7) 

Here the time evolution comes from (U), and the calcu- 
lated probability compared with experimental measure- 
ments exists only for < t < 00. 

For any experimental system present in a laboratory, 
there is also a special time: the coordinate value of the 
preparation time, tprep- It is the time on the clock at 
which the state of that system has been prepared such 
that it is representable by, say, p{t). It is also the time 



after which a detector can possibly register an observ- 
able [13: 

the observable represented by A 

is registered at f > tprep- (8) 

Comparing ([5]) with ([7]), one identifies the semigroup 
time parameter value t = of the theory with the time 
coordinate corresponding to the preparation of the state 
of an experimental system, tprep- 

Theory Physical Identification 

^ Q ^ preparation time of the state 
of the experimental system 

For non-relativistic applications, then, the relation be- 
tween coordinates and parameters is trivial. For a mea- 
surement performed at the duration t — tprep after prepa- 
ration, one has 

t t tprep- (9) 

By causality, t > tprep, and one cannot contradict what 
is in 

C. Experimental Ensembles 

The results of quantum mechanical calculations are 
probabilities meant to be compared with averages over 
many identical measurements performed on identically 
prepared, experimental systems For actual experi- 
ments, of course, one replaces the notion of "identical" 
with something like "as similar as possible." By ensem- 
ble, then, we refer to the collection of such identically 
prepared, experimental systems meant to undergo mea- 
surement. This is the usual understanding of ensembles, 
but we emphasize that we have no comment on the En- 
semble Interpretation [28| of quantum mechanics. Nor 
will we reference another level of abstraction in which 
one imagines an ensemble of state vectors f29( or of wave 
functions 3]. 

While a density operator represents in principle the 
state of a single member as well as the state of the entire 
ensemble, a measurement on just one member alone is 
useless. An ensemble always exists for any quantitatively 
useful experiment. 

Because of our identification of the parameter value 
t = with the time coordinate corresponding to the 
preparation of the state of an experimental ensemble, if 
p{t) is to represent the state of all members of an exper- 
imental ensemble, then we require the following rule: 

Rule for the application of theory 

Every member of the ensemble, at the moment it is 
prepared such that it is representable by the density 
operator p{t), is represented by that operator at the 
parametric time zero, p(t = 0), regardless of the 
coordinate time in the lab at which the preparation 
occurs- 
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This rule is already followed in the reduction of ex- 
perimental data. This is especially clear in dynamical 
experiments performed repeatedly on single members of 
an ensemble [30l - [36| . Because the rule is followed when 
reducing data, its appearance in the theory is quite nat- 
ural. 



D. Theoretical Image of an Open System 

To understand why predictions of decoherence follow 
for open systems, consider a hypothetical sequence of 
three events pictured in the timeline in Figure [TJ 

1. A physicist prepares the state of N members of 
an experimental ensemble at, according to the lab- 
oratory's clock, tprep- Depending on the experi- 
ment, of course, individual members of the ensem- 
ble may have different values of iprep- In some 
experiments, these different values are recorded as 
"time stamps." 

2. After a duration of t — tprep = 3 s, something un- 
controlled (and not directly known) perturbs some 
members of the ensemble. 

3. After the uncontrolled occurrence, the system 
evolves undisturbed until finally, at t — tprep — 8 s, 
the physicist performs an active measurement. 

Because of the uncontrolled occurrence in Figure [1] the 
system is called an open system. One can always gener- 
alize it, but for our purposes, this very simple example 
suffices. 





Evolution of a 
Quantum System 


Evolution of a 
Quantum System 













something 
preparation . , , 

uncontrolled happens 



t = Os Os 5s 

time parameter for the affected subset 



FIG. 1: Schematic timeline for an open quantum system. 

To construct the theoretical image of this sequence, 
standard quantum mechanics requires the two closed sys- 
tems pictured in Figure [2] The vectors 0„ span the space 
$0, and they can describe the state of the experimental 
system of interest. The vectors Xn span <i>^, and they can 
describe the state of the uncontrolled system. By defini- 
tion, any density operator defined in $0 can describe the 
state of the system of interest and of nothing external to 
it. 

Assume that the physicist prepares the ensemble of size 
N to be in a pure state. In the theory, one then represents 



the state of the ensemble with an operator projecting 
into a one-dimensional subspace of $0 spanned by, say 
(j)i. Before any uncontrolled experimental occurrence, 
i.e., for t — tprep < 3 s in Figure [1] the appropriate density 
operator at t is 

p{i~ tprep) ^ P^t) = U{t)\(t)i){(t)i\U\t). (10) 

In (|10p. the identification t <^ t — tprep is implicit. 

If the uncontrolled event in Figure [T] does not coher- 
ently affect the entire ensemble, it affects only a subset 
of size m. And by definition, the uncontrolled occurrence 
cannot be described by the (closed) dynamics of the sys- 
tem of interest. The state of the m affected members of 
the ensemble simply ceases to be representable by p,p(t). 
Instead, one must represent their state with a new op- 
erator, Px(t), which is in general defined in a different 
space, And according to the rule above, one must 
initially represent the state of the affected subset with 
p^{t = 0). By preparing (passively) the state anew, the 
uncontrolled event resets the time parameter, t. This is 
noted in the bottom of Figure [1] 

The state of the unaffected N — m members remains 
representable by p,f,{t). For t — tprep > 3s, 

pit -tprep) ai P^{tl) + 02 Px(*2)- (H) 

In (jlip . ai and 02 are weights that will depend on N and 
m, as well as on the nature of the physical systems. 

Note that p^{t) = Px{t) is possible but certainly not 
required. And even if they are equal, p{t — tprep) h^ (fTTj) 
will represent a mixed state because ti ^ t2 for any given 
value of the coordinate t (see Figure [TJ) 

If $0 y^ then (fTTj) may not be weU defined. To 
treat it correctly one must eventually define how an op- 
erator defined on $0 might be extended or limited to the 
space When represents an uncontrolled environ- 
ment, having $0 7^ will generally lead to predictions 
of dissipation. Further discussion can be left for another 
paper because, for the experiments we analyze in Sec- 
tion |lVl from the data we will deduce that, immediately 
after the uncontrolled occurrence, $0 = 

Now consider the time parameter for the affected sub- 
set of the ensemble. At measurement 

8a — i — tprep ^ t2 = 5s. (12) 

To conclude ([T^ . we have made use only of the defini- 
tion of open quantum systems and of the intrinsic time 



As pictured, $^ C 





set of vectors {0} = <i>^ 
Experimental system of interest 




set of vectors {x} = 
Uncontrolled system 



FIG. 2: The extent of two closed systems. 
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asymmetry; our treatment of this hypothetical sequence 
is clearly independent of any model. The time parame- 
ter values, t, do not correspond simply to the passage of 
absolute time in the physical universe. 

If the process of moving from (flUl) to (|lip can be called 
"branching" of the density operator, then after multiple 
instances of environmental interference one has a density 
operator with many branches: 

p(i- iprep) <^ ai P4,{ti) +02 Pxit2) + ... = pi{t}), (13) 

where {t} is a set of time parameters. 

A density operator for a pure state always branches 
into a density operator for a mixed state. This is a predic- 
tion of decoherence, which is an experimental signature 
of the quantum mechanical arrow of time. And though 
this branching is for the theory very simple, given our 
choice of boundary conditions it follows necessarily from 
the correct application of standard quantum mechanics. 

One could conceivably develop this approach to deco- 
herence without using the semigroup at all. Based upon 
the remarkable experiments in which single members of 
an experimental ensemble are repeatedly prepared and 
measured [SOhSO] , from physical considerations one must 
already reset the time parameter to t — whenever en- 
semble members are actively prepared anew. But we have 
included the full development here because we enjoy how 
time asymmetric boundary conditions ([3]) lead automat- 
ically to the experimental signatures of irreversibility. 

IV. APPLICATION: RABI OSCILLATIONS 
EXPERIMENTS 

As a test of our approach to open quantum systems, 
let us now analyze actual experimental systems within 
the framework of microphysical time asymmetry and the 
branching density operators. Our approach will allow 
us to address the indistinguishability of experimental en- 
semble members, and as a consequence we will find that 
our model compares very well with measurements. 

A. Rabi Oscillations 

The theory of Rabi oscillations is well established [13] ■ 
Consider a two level system, with levels described by the 
projection operators \g){g\ representing the ground state 
and |e)(e| representing the excited state. When the sys- 
tem is prepared at t = to be in the excited state, |e), 
and it is coupled to a correctly tuned radiation field, the 
Born probability at t to find the system in the ground 
state, \g), is 

V\g)ig\{p{t))^siii\nt) = i(l-cos(2m)), (14) 

where fl is called the Rabi frequency. (The probability to 
find the system in the same state in which it was prepared 



is cos^(r2t).) These expressions are calculated for closed 
systems using standard quantum mechanics. The t in 
([T4l) is therefore a single time parameter rather than a 
time coordinate or a set, {t}, of time parameters. 

In actual experiments, one never measures the result 
in p4p , because in practice, a physical system cannot be 
isolated from its environment. In this paper we wish to 
consider several very different but very clean and well- 
controlled experiments JTI, [HM llUSi ■ 

In one experiment [34l], internal levels of a ^Be"*' ion 
couple to a harmonic binding potential. Rabi oscillations 
occur between two of the coupled internal and vibrational 
levels, and the oscillations between different sets of lev- 
els are measured. In another experiment [s^, Rabi os- 
cillations are observed between the circular states of a 
Rydberg atom coupled to a field stored in a high Q cav- 
ity. In 36], the Rabi oscillations are between the spin 
states of two electrons in a double quantum dot. In [38| . 
Rabi oscillations occur between motional states of elec- 
trons bound to shallow donors in semiconductors. These 
motional states mimic the levels of a hydrogen atom. In 
three experiments (l^ . [sol . |40| , oscillations occur between 
excitation levels of electrons confined in quantum dots. 

In [33 - [36| . one does not match the prediction in (|14l) . 
Instead, the measured probability is fit by an appropri- 
ately damped sinusoid, 

^|9>(9|(P(*)) = ^(1 - e-^*cos(2r!?)), (15) 

where 7 is an experimentally determined damping factor. 
Dynamical measurements are made in the lab, and the t 
in (jlSp is therefore a time coordinate. For simplicity, we 
have begun to set tprep = 0, with the understanding that 
t hereafter represents a duration from tp^ep- 

In [36} , amplitude, offset, and phase are fit as well, 
because the oscillations are measured also as a function 
of a swept detuning voltage. In [ll, [i^, photocurrent 
is measured, and from it the decoherence of the Rabi 
oscillations is inferred. Numerically solving the models 
they fit to their data, the subsystems undergoing Rabi 
oscillations fit ([15]). 

Experiments also reveal that the damping factor, 7, 
depends generally on the Rabi frequency, fi, regardless 
of the nature of the experiment. In the experiment with 
^Be^ ions [34], the different levels are described by the 
kets \ i,n) and ] f, + 1), where ] J,) and ] f) are internal 
states of the Be ion, and \n) represents vibrational Fock 
states. Rabi oscillations are measured for the frequen- 
cies [13, Ell 

r!„,„+i - n- ; (0.2022), (16) 

V" + 1 

where L\ is the generalized Laguerre polynomial. The 
corresponding damping factor, 7„, is measured to in- 
crease with n according to 

^ w (1 nf-\ (Measured) (17) 
70 
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In the experiment with spin states of two electrons 
in a double quantum dot, the damping factor, 7, is stated 
to be proportional to the Rabi frequency, ft, though no 
mathematical relation is given. 

In IZ, .38r.40j , this dependence has been called Exci- 
tation Induced Dephasing (EID), and its cause has re- 
mained an open question. In l40j . the dependence 
of the damping factor on the Rabi frequency has been 
found to be 7 cx ri^. 

These results for such different experiments strongly 
suggest that some dependence of 7 on O is general and 
independent of experimental specifics. 



B. Development of a Model 

Our task is to develop a single model of decoherence 
that reproduces ((T5|) for systems undergoing Rabi oscil- 
lations, and, if possible, the general dependence of 7 on 
il, or EID, also mentioned above. Because of the distinc- 
tion between parameters, t, and coordinates, t, measured 
results made in coordinate time must be compared with 
probabilities calculated also in coordinate time. We will 
therefore call VA{p{t)), which is a function of coordinate 
time, the "predictive probability." And below we will 
take A = \g){g\. 

Using the branched density operator from (|13p . we 
have the predictive probability 

VA{p{i)) = Tr(Ap(i))^Tr(Ap({t})) 

= Tr(Aaip^(ii)+Aa2/3x(*2) + ...) 

= airA{p4ti)) + a2rA{pxit2)) + ■ ■ ■ (18) 

The long-time behavior of the measured result in (jlSp 

'^\a){a\{pi^)) = ^|e>{e|(p(*)) \- From this we can 
deduce that there is no dissipation for the experimen- 
tal system, and that the term VA{Px{i2)) in ought 
not to introduce losses. For these clean experiments, we 
then hypothesize that immediately after the uncontrolled 
occurrence, $0 = Therefore, we want a predictive 



probability of the form 



VA{prn{tm)) (19) 



where Mhranch is the number of branches. Note that 
we no longer require the label for p. This formula is 
deceptively simple in appearance, because, in general, the 
number of branches, Mtranch, will depend on t, and the 
value of am' for a fixed m' will depend on Mbranch- 

The time parameter values tm are related unambigu- 
ously to t via the rule in Section flll CI The a™ will come 
from models. 

We still require a prescription for the pm- Leggett has 
explained [i^l that decoherence "is exactly the result of a 
'measurement' whose result is uninspected," and he has 
called this process "garbling." Measuring a system to be 
in a state is equivalent to preparing that system to be 



in that state, so these passive measurements that garble 
an ensemble are also passive preparations. We will there- 
fore assume that the effect of an uncontrolled interference 
event is to measure and therefore prepare anew some sub- 
set of our ensemble. And for obvious reasons, one should 
assume that nature does not distinguish between passive 
and active preparations. 

A hypothesis physically equivalent to garbling is that 
the uncontrolled environment interacts only with energy 
eigenstates of the system of interest. Then, because = 
$0, when the uncontrolled event occurs at tm = 0, the 
most general form for pm is 



Pm{tm = 0) = w„i \g){9\ + Vm |e)(e|. 



(20) 



with Wrn + = 1. In (120^ . w is a weight appropriate 
for ensemble members prepared passively to be in the 
ground state, and w is a weight for ensemble members 
prepared in the excited state. The weights are of course 
independent of tm, and their values will also come from 
our model. 



C. Predictive Probability 

Given our hypotheses and deductions, the simplest rea- 
sonable and general model the author can conceive is: 

1. Members of the experimental ensemble are actively 
prepared in the lab at t = tprep- Again, for simplic- 



ity, we set tp 



0. 



2. The members of the ensemble can possibly suffer 
environmental interactions at the times nAt, where 
n — 1,2,3... As a result, some systems are pas- 
sively prepared to be in either the ground state or 
the excited state. We will call these interactions 
"interference events." 

3. At every nAt, there is some probability, (1 — 77), for 
a member to suffer perturbation and thereby to be 
prepared passively. 

In the second step, we have introduced the time scale. 
At, of interaction with the environment. 

In the third step we have introduced a parameter, 77, 
with < 77 < 1, to represent the susceptibility of en- 
semble members to environmental interference. As de- 
scribed, the parameter rj is the probability for any given 
member of the ensemble not to suffer interference at one 
of the times nAt. For a perfectly isolated system, 77 = 1. 

According to this simple model, we can write a general 
formula for the predictive probability as more and more 
branches form: 



^19>(91(PW) = ' 



Po{i) 0<i<lAi 
pi{i) lAi<i<2Ai 



(21) 



Pn{i) nAt <i < {n + l)At 
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Here, every Pn{t) will have the form of the sum in (flQl) . 
and Mbranch = n. 

Let us assume that a physicist identically prepares ev- 
ery member of the experimental ensemble initially to be 
in the excited state, |e). Because no members will have 
suffered environmental interference before t — 1 At, we 
have for the initial value po(i) = sin^(f2f). 

At t = 1 At, the fraction (1 — r/) of the ensemble will 
be passively measured and thus prepared in a new state, 
with new initial conditions. At this point, Mi,j.anch — Ij 
and from (IT51) we have 

Piii) ^ aoV\g^g\{po{to)) + aiVig^g\{pi{t,)). (22) 

The first term in (1221) . oq 7^|g)(g| (/9o(^o)) i corresponds to 
the subset of the ensemble that passes the time lAt with- 
out suffering interaction with the environment. Accord- 

I 



ing to our model, it is r]sm^{ni) — r]po{i). For this 
unperturbed subset, of course, t <^ tg. 

The second term in (f22|) . aiV\g)(g\(pi{ti)) , is the first 
branch. It corresponds to the subset of the ensemble that 
suffers perturbation and is thereby garbled, or passively 
prepared anew. According to our model, ai = (1 — 77). 

By our prescription for the density operators, we know 
that pi{ti = 0) is fixed by the state of the perturbed 
subset as it is prepared a.t i = lAi by the uncon- 
trolled event. The perturbed subset is passively mea- 
sured at t — lAt, so assuming that garbling is a fair 
measurement, pi{ti = 0) will have the form of (I^Hl) with 
wi = sin^(l nAi) = po{lAi) and vi = 1 - po{lAi). 

Finally, with the identification ti t — lAt for the 
perturbed subset (see Figure [ij, it is straightforward to 
find from the theory of Rabi oscillations 



Pi 



(i) = ripoii) + il-ri)(po{lAi)cos^{n (i - lAt)) + (l -po(lAi)) sm^{n {i - lAi)) 



(23) 



r 



This is a nice form because we can stop writing out probability. For general n, we have 
branches and use a recursive formula for the predictive 

I 



Pn{i) = VPn^i) + {1 - T])(pn-i{nAi) cos^ {n (i-nAi)) + (l - p„_i(nAt)) siu^in (<- nAi))). 



(24) 



Figure [3] shows the results of two sample calculations 
using ([21) with (|24|) . The dots are from the recursive 
calculation using our model. The solid lines are plots of 
the decaying sinusoid in (jlSp . which fits the experimental 
measurements. For both figures we have used fiAt « 
0.08. The results in Figure 3(a) were calculated using 
77 = 0.99 and resulted in a fitted value for the damping 
factor of -f/n = 0.05. The results in Figure [3(b)] were 
calculated using ry = 0.997 and resulted in a fitted value 
for the damping factor of = 0.015. Recall that rj — I 
for a perfectly isolated system. 

Our calculation of the predictive probability results in 
no frequency shift away from the Rabi frequency, fl, at 
early times, in agreement with experiments. (With the 
master equation, one predicts a frequency shift and must 
assume very strong driving to get rid of it. See Sec- 
tion IIVEI ) Note, however, that our method results in 
a fitted damping factor, 7, that is independent of ft. 
Though our first attempt has been instructive, it is thus 
not correct. We have included this incorrect result not 
only because it clearly illustrates our new method, but 
also because the recursive structure of (j24[) will reappear 



in the correct model below. 



D. Predictive Probability and Indistinguishability 



Consider the experiments [3J-|36j in which one repeat- 
edly prepares and then measures single members of the 
experimental ensemble. One assumes that performing 
a number of consecutive measurements, each done on a 
single member, is equivalent to performing a measure- 
ment on the same number of non-interacting, simulta- 
neously present members of an identically prepared en- 
semble. This is a statement of ergodicity. And when the 
identically prepared members are present simultaneously, 
one treats them as indistinguishable. In our first model, 
we did not treat them so. And to enforce ergodicity, we 
would like to. 

Though in these experiments, one can easily label any 
member of the ensemble with the value of the time coor- 
dinate, t, at which it was actively prepared or measured, 
coordinate values are always physically meaningless. The 
only physically significant label for identically prepared 
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5 10 15 20 25 30 
Dimensionless time coordinate, fl? 

(a)With J? = 0.99, we fit 7/Q = 0.05. 




5 10 15 20 25 30 
Dimensionless time coordinate, fit 

(b)With -q = 0.997, we fit 7/C = 0.015. 

FIG. 3: The predictive probability, "PjgXgi when mem- 

bers of the ensemble are distinguishable. For both plots, 
i}At « 0.08. The dots are the results of recursive calcula- 
tions using (|24|) and (|21|l . The solid lines are plots of the 
damped sinusoid (|15|) that fits the experimental data. Recall 
that rj is the probability that members will not suffer a per- 
turbation at the times nAt. Values of rj and At were chosen 
for aesthetics only. These plots do look good, but our results 
are not yet correct. See the text. 

members is the duration between preparation and mea- 
surement. And, even in these experiments, one cannot 
in principle label any given member of an ensemble with 
a duration. To understand why, consider the result of a 
hypothetical, active measurement on any single member 
of an experimental ensemble: 

1. At t — a physicist actively prepares the members. 

2. At i — lAi subset A is passively prepared to be in 
a mixture of the ground state and the excited state. 

3. At t — 2 At subset B is passively prepared to be in 
a mixture of the ground state and the excited state. 

4. At some time coordinate value t > 2At a physi- 
cist actively measures a single member to be in the 



ground state, \g). 

Because of the probabilistic nature of quantum mechan- 
ics, it is impossible to distinguish if the actively measured 
member was in subset A, subset B, both, or neither. 

In other words, if the state of an ensemble is repre- 
sented by the branched density operator 

p{{t}) = ao paito) + UA PA{tA) + as PBits), (25) 

then it is impossible to know to which branch any given 
member belongs. 

One way to distribute indistinguishable ensemble 
members among several branches is to count different 
combinations of the members themselves. But we can- 
not use such an approach for two important reasons. The 
first is because quantum mechanics makes predictions for 
ensemble averages, and density operators simply do not 
follow the trajectories of individual members. The sec- 
ond reason is that any model based on such counting will 
depend on the size of the experimental ensemble. This is 
clearly unacceptable. 

If, however, one imagines that ensemble members are 
labeled by the individual interference events at which 
they have been most recently prepared, we hypothesize 
that we can enforce ergodicity by treating the events 
themselves as indistinguishable. The resulting models 
will depend not on the ensemble size but instead on the 
number of interference events, which is exactly what we 
want. 

When modeling experiments, one would probably use 
a full simulation. To maintain transparency, however, we 
will again search for an analytical formula. First, assume 
that we wish to calculate the predictive probability given 
that there have been n chances for ensemble members 
to have suffered a single interference event. If interfer- 
ence events can possibly occur at intervals of At, the pre- 
dictive probability can be written V\g)(^g\(^p{nAt)'j . Our 
scenario is as follows: 

1. Choose a combination of time intervals (of length 
At) such that systems have survived only a dura- 
tion of 1 At before being passively measured. 

2. Adjust the probability at nAt accordingly. 

3. Because the systems and interference events are in- 
distinguishable, so are the intervals between the 
events. Put the intervals chosen above back into 
the original set. 

4. Choose a combination of time intervals such that 
systems have survived a duration of 2 At before be- 
ing passively measured. 

5. Adjust the probability at nAt accordingly. 

6. Put the intervals chosen above back into the origi- 
nal set. 

7. Repeat until reaching nAt. 
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We have again chosen to parametrize our model with 
a time scale, At. However, we will require a slightly 
different interpretation for the number parameterizing 
the members' susceptibility to environmental interfer- 
ence. We shall now use /3 (rather than 77), with < /3 < 1 
and (3=1 for a perfectly isolated system. Physically, (3 
will be the probability that a randomly chosen time in- 
terval will have come before an interference event. This 
can be better understood after equation ((3T|) . 

To implement steps [T] and HI above, we will make use 
of the binomial distribution, 

5(n,fc,/3)^ Q/3'=(l-/3)"-^ (26) 

The distribution in (|26p gives the probability for the oc- 
currence of any combination, regardless of order, of k 
events with probability /3 and (n — k) events with prob- 
ability (1 — /3). Because the distribution gives probabil- 
ities for combinations rather than permutations, with it 
we can treat indistinguishable interference events. Note 



The explanation of is straightforward. We need 
to relate the binomial distribution to the passage of time, 
so at every step k, with < k < n = -i, we will count the 
(normalized) number of combinations for arranging the n 
time intervals such that k of them came before the single 
interference event. This number is given by b(ri,k,/3). 
Because possible interference events occur at increments 
of the time scale. At, the weight b{n, k, [3) must then be 
attached to any passive preparation occurring at kAt. 

Again, branches of a density operator form when sub- 
sets of the ensemble suffer environmental interference and 
are passively prepared. For a branched density operator, 
in each line of ([^5]). the term multiplied by b{n,k,(3) is 



also the normalization 

E(^)/3'=(1-/3)"-' = 1. (27) 

Using ([27)1 we will relate the binomial distribution to a 
probability. 

To implement steps [2] and [5] in this model, we will 
once more use the rule deduced in Section UlI C[ to reset 
to < = the time evolution parameter of the density 
operator corresponding to members passively prepared 
as a result of environmental interference. 

To see how the binomial weights work, let us write 
a simple formula for the case that individual ensemble 
members will have suffered only one interference event 
before an active measurement occurs. This will truncate 
our formula at a reasonable size, and we will generalize 
to multiple events below. The probability at AAi to find 
in |g) members that have been initially prepared at t = 
to be in |e) is 



+ 
+ 
+ 
+ 

(28) 



the Born probability for the branches formed at kAt. 

Consider the second line, where fc = 3, and compare 
to the notation in (gHl)- The weight w = sin^(fi3Af). 
The weight v = cos^(ri3Ai). The duration from passive 
preparation to 4At is then (4 — 3) At = lAi, so the func- 
tions with argument lAi contain the time dependence of 
the Born probability of the mixed state. 

Let us introduce the notation V^g^i^g^[p{nAi)) to rep- 
resent the predictive probability under the assumption 
that members on average will have suffered i or fewer in- 
terference events before measurement. Then for general 
n, 



I 

P|g)(g|(p(4At)) = 6(4,4,/3)(sin2(r2 4Ai)cos2(17 0At) -|-cos2(r2 4Ai)sin2(r2 0Ai) 

6(4, 3, (3) (^sin^ifl 3Ai) cos^{n lAt) + cos^{fl 3 At) sin2(r2 lAi) 

6(4, 2, (3) (sin2(rj 2 At) cos^{n 2 At) + cos^{n 2 At) sin^in 2 At) 

6(4, 1, (3) (sin2(r2 lAi) cos^{n 3 At) + cos^{n lAi) sm^{n 3 At) 

6(4, 0, 13) (^sin2(r2 OAt) cos2(f7 4 At) -I- cos^{fl OAt) sin2(r2 4At) 
I 



I 

n 

'^\9){9\ (P'y^^^')) = /3)(sin^(^ kAi) cos^{n{n ~ k)Ai) + cos^{n kAt) &m^(n(n - k)Ai)^ . (29) 

k=0 

I 

By simply exchanging the cos^(f2fcAt) and sin^(f7fcAt) terms, we calculate T^.^V, , (p(nAt)) , which is the proba- 
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bility to find the ensemble members in the excited state. 

In we have assumed that systems will have suf- 
fered at most one interference event. To allow for the 
possibility of multiple events, the terms sin^(f2/cAi) and 



cos^(ri kAt) must be replaced with new functions of kAt, 
that predict the effects of interference events prior to the 
single event assumed in (f29l) . For general i, 



n 



fc=0 



In (PU)) . the recursive structure of (dJ) and ([^ has 
reappeared, but here the ensemble members and inter- 
ference events are indistinguishable. Handling the possi- 
bility for more and more interference events requires the 
nesting of more and more summed terms into pop . These 
nested equations are significantly more difficult to solve 
than are the recursive equations used for distinguishable 
events. 

We have written the predictive probability as a func- 
tion of nAt. The final step is to scale our result back to 
the time coordinate, t. The first moment of the binomial 
distribution is 



(31) 



k=Q 



After stepping through time to nAt, on average /? n of the 
intervals will have preceded the interference event num- 
ber i. This provides us a physical interpretation of our 
two parameters, and to ensure that the time scales with 
something physical, we will need to use {k)Ai = (3 nAi = 
t. After a calculation of the predictive probability as a 
function of nAt, we make the replacement 



t 



pAt 



(32) 



This restricts us to non-zero values of (3 and At. We have 
also simply interpolated between the discrete values of 
n at which (|30p is actually defined. Because our time 
scale is understood to be an average value, it would be 
inappropriate to assume for our model anything more 
complicated. 

The predictive probability at t, assuming i interference 
events, to find in the ground state a system initially pre- 
pared in the excited state is therefore 



(^) 



t 



13 At 



(33) 



Numerical solution of (1501) is straightforward. One can 
also perform the summations and find closed form expres- 
sions for the predictive probabilities. In Figure 2] we have 

,(5) 



fewer interference events. For clean experiments, with (3 
close to 1, we expect a good approximation for i = 5. 
Again, the agreement with the experimentally measured 
damped sinusoid is very good. 




5 10 15 20 25 30 
Dimensionless time coordinate, £Lt 

FIG. 4: Plot of the predictive probability, 7'|3)(gi (p(i)) , for 
indistinguishable ensemble members. For the dots we have 
used equations 1)30^ and (|33p . The solid line is a plot of the 
damped sinusoid (|15|l that fits the experimental data. We 
have used SlAt ~ 0.7 and /3 = 0.995, and we have fit 7/^ = 
0.039. 

As opposed to the previous results for distinguishable 
ensemble members, however, when we treat the interfer- 
ence events as indistinguishable, we find that the damp- 
ing factor, 7, does in general depend on the Rabi fre- 
quency, f2, in agreement with measurements and the phe- 
nomenon called Excitation Induced Dephasing. Further- 
more, in our simple model we have assumed only a time 
scale. At, and a parameter, /?, describing the experimen- 
tal system's susceptibility to interference from the envi- 
ronment. It can be applied to any real system undergoing 
Rabi oscillations. 

For j3 close to 1 , the predictive probability is dominated 
by the terms proportional to bin, k = n, /3). Solving the 
rather crude, truncated form 



plotted r^g^ (pit)) , which is the predictive probability ^^^^ (p(nAf)) « ^ 6(n, fc, /3)sin2 {Q kAt) (34) 



given that there have been for the average member 5 or 
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and using (l33l) . we get 



-(2-(l-/3(l-e 



))" 



-(l-/3(l-e+2^^*"))^).(35) 



Finally, we can conclude that a dependence of 7 on fJ, 
or BID, is indeed general, as suggested by experiment, 
and that it is a measurable effect of the indistinguishabil- 
ity of separate, uncontrolled interactions between quan- 
tum systems and their environment. 



For small At, (|3S]) is 



E. Comparison With the Standard Program 



,-7t 



(cos(2f7f) + 0(AP))), (36) 



where 



7 = 2(1-/3)172 a?+0(AF 



(37) 



For /3 « 1 and ilAt <C 1,7 is quadratic in ft. This 
matches the results in [H, [io] ■ And when the systems 
are treated as indistinguishable, we require none of the 
typical assumptions reg arding the frequency response of 
perturbation terms [l2l l40l lisl 13] ■ 

The same model is general enough to apply also to the 
system in 3J], and with it wc fit the measured relation 
p?)) between 7 and SI. Figure [5] is the result of fitting 



the damping factor, 7„, to 'P\g){g\{pi^)) ^ calculated with 
the sequence of frequencies in (fT6)) and with iloAt w 0.2. 
(The exponent can be shifted by choosing different time 
scales.) Unfortunately, we are limited by computational 
resources to n < 6. 



O 

■4— » 

B 
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FIG. 5: Matching the experimental results for the ratio of 
damping factors, The large dots result from our theoret- 
ical calculation of the predictive probability using H30p and 
([33}. The solid line is the experimentally measured relation, 
(1 + n)'^''^ . To indicate a scale for the exponent, the thin 
dashed line is a plot of (1 -I- n)" *, and the thin dotted line is 
a plot of (1 + 71)" ^ 



Figure [5] shows a very good quantitative agreement 
with experiment (jl7p , and again we have not required any 
experimentally specific assumptions. Our study indicates 
that the decoherence measured in [3^ results from envi- 
ronmental interference having a characteristic frequency 
of sa 5Slo- 



In the standard decoherence program, one treats these 
experiments using the decoherence master equation. It 
is the method of choice when one works without time 
asymmetric boundary conditions and when one identi- 
fies t with coordinate time. The generic solution of the 
master equation describing a system undergoing Rabi os- 
cillations and on resonance is 



jAIE 



4n^ 



r2 + 81^2 



e"^^*/^(cos fit + ^sin fit)) . 



(38) 

To match our convention, we have added a factor of 
2 to the definition of Rabi Frequency in In ([38]), 

fi = '\J 4172 — (j)^. Unless a specific environmental inter- 
action has been assumed, F is the spontaneous emission 
rate for the system in its excited state. It is a constant, 
independent of the Rabi frequency. 

A significant shift from the Rabi frequency, 17, is not 
observed in experiments, however, so one must assume 
very strong driving: 2S7 ^ r/4. The limit of ([55]) is then 



3-3rt/4cos2m) 



(39) 



One can see immediately why the measured results 
have been puzzling. For general calculations, the solution 
_of the master equation predicts for the experiment 

m 



70 



1, 



(Master equation) 



(40) 



and that 7 is independent of Vt. 

In the hope of matching experiment, one typically as- 
sumes a detailed form for the interaction operators. For 
the experiment in [s^l , in which (|17l) was measured, sev- 
eral such studies have been carried out ;45-5ll. None of 
these studies, however, has resulted in an agreement as 
quantitatively good as that in Figure [5j 

Furthermore, the models for decoherence in [45l - [5l| 
have necessarily been highly tuned. They are thus not 
applicable to the other types of experiments [l^ |H, Issl - 
lioj . And, with the exception perhaps of [13], they sug- 
gest that a dependence of 7 on SI is not general, in 
disagreement with measurements. Similarly, the mod- 
els [13, SO, in, nil used to explain the EID measured in 
experiments with shallow donors and quantum dots can- 
not be applied to ions in a Paul trap 34] . 

This paper describes the general framework that one 
must use when applying a time asymmetric theory. Full 
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treatment of experiments will likely require more de- 
tailed models addressing multiple sources of environmen- 
tal interference, dressed states, etc. But the qualitative 
and quantitative success of our single, general model for 
different types of Rabi oscillations experiments is very 
promising. 

V. CONCLUSION 

A choice of boundary conditions results in intrinsic 
time asymmetry endowed at the microphysical level, even 
for closed systems. The theoretical expression of this 
asymmetry is time evolution generated by a semigroup. 
Constructing the theoretical image of open systems re- 
quires a new understanding of how time parameters cor- 
respond to the passage of time in the physical universe. 
As a result, standard quantum mechanics already pre- 
dicts decoherence, without invoking a reduced dynamics. 
This suggests that the extrinsic arrow of time may be 
the experimental signature of a more fundamental, mi- 



crophysical time asymmetry. 

The practical result is a new framework for the treat- 
ment of decoherence, in which states of open systems are 
represented by branching density operators rather than 
by solutions of the decoherence master equation. As an 
application, we have created a simple model matching 
experimental results from Rabi oscillations experiments. 
With the new formalism, we can conclude that a general 
yet puzzling experimental result, known as Excitation In- 
duced Dephasing, is the measurable consequence of the 
indistinguishability of separate, uncontrolled interactions 
between quantum systems and their environment. 

The new formalism is very promising for the study of 
quantiim decoherence. In forthcoming papers, we will 
show how the framework naturally extends to statistical 
mechanics and the increase of quantum mechanical en- 
tropy. And perhaps even more compelling, we will also 
demonstrate that this same framework that works so well 
for decoherence is also very effective when applied to scat- 
tering theory. 
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